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ABSTRACT 

High-resolution X-ray observations have revealed cavities and 'cold fronts' with sharp edges in tempera- 
ture and density within galaxy clusters. Their presence poses a puzzle since these features are not expected 
to be hydrodynamically stable, or to remain sharp in the presence of diffusion. However, a moving core or 
bubble in even a very weakly magnetized plasma necessarily sweeps up enough magnetic field to build up a 
dynamically important sheath; the layer's strength is set by a competition between 'plowing up' and slipping 
around of field lines, and depends primarily on the ram pressure seen by the moving object. In this inherently 
three dimensional problem, our analytic arguments and numerical experiments show that this layer modifies 
the dynamics of a plunging core, greatly modifying the hydrodynamic instabilities and mixing, changing the 
geometry of stripped material, and slowing the core through magnetic tension. We derive an expression for the 
maximum magnetic field strength and thickness of the layer, as well as for the opening angle of the magnetic 
wake. The morphology of the magnetic draping layer implies the suppression of thermal conduction across the 
layer, thus conserving strong temperature gradients. The intermittent amplification of the magnetic field as well 
as the injection of MHD turbulence in the wake of the core is identified to be due to vorticity generation within 
the magnetic draping layer. These results have important consequences for understanding the complex gasdy- 
namical processes of the intra-cluster medium, and apply quite generally to motions through other magnetized 
environments, e.g. the ISM. 

Subject headings: hydrodynamics — magnetic fields — MHD — turbulence — galaxies: clusters: general — 
diffusion 
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1. INTRODUCTION 

Recent observations of very sharp 'cold fronts' in galaxy 
clusters raise unanswered questions in the hydrodynamics of 
galaxy clusters (see for instance the review of Markevitch and 
Vikhlinin 2007), for such abrupt transitions are not expected 
to be stable against either hydrodynamical motions or diffu- 
sion for extended periods of time. 

It has been known for some decades in the space science 
community that an object moving super- Alfvenically in a 
magnetized medium can very rapidly sweep up a significant 
magnetic layer which is then 'draped' over the projectile {e.g., 
Bernikov and Semenov 1980). For concreteness in discussing 
the process, we show in Fig. 1 a picture of this mechanism 
taken from one of our simulations, which will be described in 
more detail in later sections. 

There has been significant interest in applying this idea 
of magnetic draping in galaxy clusters {e.g., Vikhlinin et al. 
2001; Lyutikov 2006; Asai et al. 2004, 2005, 2007) as such 
a magnetic field could naturally inhibit thermal conduction 
across a front {e.g., Ettori and Fabian 2000) allowing it to 
remain sharp over dynamically long times. Although such 
draping has been explored in the past, in the space sciences 
the resulting dynamics is relatively simpler, as generally the 
object being draped is a solid body, with little interior dynam- 
ics of its own. However, in the case of for instance a merger 
of gas-rich clusters, the hydrodynamics of the draped plung- 
ing core can also be modified, with the strong magnetic layer 
providing some stabilization against instabilities that would 
otherwise occur (Dursi 2007). 

The effect of a strong draped magnetic layer could be even 
greater for underdense objects, such as for bubbles moving 
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through the intercluster medium, as seen at the centers of 
many cool-core clusters {e.g., McNamara et al. 2005; Birzan 
et al. 2004). In this case, the bubble would be quickly dis- 
rupted on rising absent some sort of support {e.g., Robin- 
son et al. 2004). However, the draping of a pre-existing 
magnetic field may strongly alter the dynamics and suppress 
hydrodynamic instabilities, as seen recently in simulations 
(Ruszkowski et al. 2007b). The morphology of the draped 
magnetic field may be able to suppress transport processes 
across the bubble interface such as cosmic ray diffusion and 
heat conduction. This has important consequences for cosmic 
ray confinement in these buoyantly rising bubbles and may 
explain re-energized radio 'relic' sources, broad central abun- 
dance profiles of clusters, and the excitation of the Ha line in 
filaments trailing behind bubbles (Ruszkowski et al. 2007a). 
Although the analytics and simulations we discuss here fo- 
cus on the case of an overdense 'core' moving in an external 
field, we expect the basic magnetic dynamics to also extend to 
the case of an underdense bubble probably depending on the 
magnetic energy density of the plasma. 

Because we are interested here in the fundamentals of a ba- 
sic process — that of the draping of a field around an object 
and the resulting hydrodynamical effect on the object and its 
interaction with the external medium — we consider for this 
paper, in both our analytic and computational work, the sim- 
plest possible case; an overdense, non-self-gravitating 'blob' 
moving through a quiescent medium with a magnetic field 
uniform on the scales considered. (The term we will use for 
this blob will depend on the situation; when discussing astro- 
physical implications, we will speak of 'cores' or 'bullets', 
depending on the circumstances, or 'bubbles' for underdense 
regions; for the case of our numerical simulations, we will re- 
fer to 'projectiles', as the overdense fluid in the simulations 
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differs in structure from cores or bullets in lacking self grav- 
ity; in our analytic work where the blob is a rigid sphere, we 
will refer to the blob as a sphere or spherical body.) 

We further consider the case of the object moving subson- 
ically; while the case of supersonic motion is interesting and 
highly relevant, we anticipate that in the usual case where 
the bow shock is well separated from the magnetic layer - 
that is where the standoff distance A w '^R ^ R, (where R 
is the radius of the core, po is the ambient density, and ps is 
the shocked density) is much greater than the magnetic layer 
thickness I « ■M'a^R ^ R, (where I is the approximate mag- 
netized layer thickness and Ma the alfvenic Mach number, 
as discussed in more detail in the next section) that the argu- 
ments here will also hold, so we save the more complicated 
geometry and larger parameter space of the compressible case 
for future work. 




Fig. 1 . — A rendering of one of the three-dimensional simulations (referred 
to as Run F later in this work) performed for this work, discussed in more de- 
tail in later sections but included here to illustrate the physical picture. An 
overdense projectile is sent through a uniformly magnetized medium, sweep- 
ing up magnetic field ahead of it. Plotted is a density isosurf'ace, coiTespond- 
ing to the mean density of the bullet, and some fiducial magnetic field lines. 
The cut-plane is coloured by magnetic energy density, as are the field lines. 
The minimum of the colourmap for the magnetic energy density is set to be 
the magnetic energy density in the ambient medium; energies less than that 
(which occur in the wake, due to entrainment of unmagnetized material as the 
bubble rises) is set to the color of the ambient energy density. The magnetic 
field is 'draped' into a thin layer forming a bow wave, leaving turbulence in a 
wake behind the bullet. Magnetic field lines pile up along the stagnation line 
of this initially axisymmetric bullet, while in the plane perpendicular to the 
initial field, the field lines can slip around the bullet. A PDF version of this 
manuscript with an interactive 3d version of this figure, following Barnes and 
Fluke (2007), is available at http://www.cita.utoronto.cay~ljdursi/draping/. 

In §2 we give an overview of the physics of draping, putting 
our work in the context of previous results; in §3 we describe 
our analytic and computational approaches; in §4 we compare 
the results of our two approaches, and from the understanding 
gained there in §5 we describe characteristics of draping; we 
discuss the effect of instabilities in §6, consider the limitations 
of our results and consider applying them to later times in §7, 
and finally conclude in §8. 

2. GENERAL PHYSICS OF MAGNETIC DRAPING 

Previous work {e.g., Bernikov and Semenov 1980; Lyutikov 
2006) has looked at the basic picture of magnetic draping 
in a simplified way in some detail; we summarize some of 
their key results as well as our new insight into this problem 
here. In these works, the known potential flow around a solid 
sphere is taken as an input, and a purely kinematic magnetic 
field, uniform and perpendicular to the direction of motion, 
is added. The derivation of Bernikov and Semenov (1980) 
is clarified, and a novel set of useful approximations for the 
resulting field near the solid sphere are given, in Appendix A. 



Because in this case the flow falls quickly to zero at the 
surface of the moving sphere, magnetic field rapidly 'builds 
up' around the projectile, and in the kinematic limit eventu- 
ally becomes infinite. The high degree of symmetry along the 
stagnation line (the axis of symmetry of the object pointing 
in the direction of motion) greatly simplifies the mathematics, 
and as shown in for instance Lyutikov (2006), the magnetic 
field strength directly along the stagnation line is given by 

where Bo is the ambient magnetic field, po is the ambient den- 
sity, R is the radius of the solid sphere projectile, and s is 
the distance along the stagnation line from the surface of the 
sphere. 

The analytic works cited, and presented here, considered 
purely incompressible flow; for our simulations, we consider 
only very modest compressibility, with projectile motions 
through the ambient fluid quite subsonic, so it suffices for the 
moment to consider in the external medium p= po- The fluid 
here is further considered to be infinitely conducting; how- 
ever, the buildup of magnetic field without a corresponding 
buildup of mass does not violate the 'flux-freezing' condition, 
as shown in the cartoon Fig. 2 as incoming fluid elements are 
'squished' along the sides of the incoming sphere, so that the 
magnetic flux coming out the sides of the fluid element re- 
mains constant, even as the concentration of field lines builds 
up along the stagnation line. Further increase in magnetic en- 
ergy comes from the stretching of field lines in the direction 
of motion of the core. 




Fig. 2. — A cartoon showing the distortion of incoming fluid elements and 
stretching of field lines as a red spherical projectile moves upwards through 
the ambient medium. 

In reality, of course, the magnetic field does back-react onto 
the flow, and the kinematic potential flow solution fails for 
two reasons - buildup of a strong magnetic field layer (which 
violates the kinematic assumption) and creation of vorticity 
(in conflict with the potential flow assumptions). 

The magnetic field should exert a significant back-reaction 
when the resulting magnetic pressure is comparable to the ram 
pressure of the incoming material: B^/87r ^ pou^, where u is 
the speed of the core through the quiescent ambient fluid. The 
first place this will happen is along the stagnation line, which 
by symmetry will be the location of the largest magnetic en- 
ergy density. The layer of magnetic field with this magnitude 
is expected (from Eq. 1 and assuming / <C /?) to be of thick- 
ness 

l = -r^'^ (2) 

where Ma = u /va is the Alfvenic Mach number of the core, 
v\ = Bl/4-TTpo is the ambient Alfven speed, and a is the con- 
stant of proportionality describing the maximum magnetic 
pressure in units of the incoming ram pressure, B^^^^/d^TT = 
ap()U^. We will see that a « 2 and fiducial values for the situ- 
ations considered here will involve Aij^ w 3, so that a typical 
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Fig. 3. — A cartoon showing the expected geometry of the draped mag- 
netic field (blue) over the object (red). Seen in the plane of the direction of 
the ambient field, panel (a), with the direction of the ambient field shown, a 
clean bow wave is presented with a well-defined opening angle. In the plane 
perpendicular to the ambient field, panel (b), the field lines can slip around 
the projectile, and the flow would close back in on a stagnation line on the 
other side of the object except for largely-2d vortical motions induced by in- 
stabilities at the magnetic interface. The geometry of the flow in the region 
indicated by dashed box depends heavily on the final shape, and thus internal 
structure, of the moving object. 



value for / will be approximately /?/36. Even such a very thin 
layer can have important effects, both in terms of suppress- 
ing thermal conduction (Ettori and Fabian 2000) and hydro- 
dynamic instabilities (Dursi 2007). 

It should be noted here that when we use the Alfvenic Mach 
number M.a through this work it should really be considered 
a dimensionless ratio of ram pressure to magnetic pressure 
{M.\ = pQir /{2Pb„)), or at the least, some caution should be 
used when interpreting it as a ratio of velocities {u/va) as the 
velocities are oriented in different directions; in the work pre- 
sented here, the velocity of the draped object will always be 
completely orthogonal to the ambient direction of propagation 
of Alfven waves. Thus there is an important sense in which 
our projectiles are always (infinitely) super-Alfvenic, which 
is not captured in the ratio Ma- 

Sweeping up such a magnetic field will occur on a timescale 
t/tc ^ ^{1/R)Ma ~ {s/aMAT^, where tc = 2R/u is the pro- 
jectile's own crossing time. This result means that, because 
the magnetic layer is very thin, a strong field can be built 
up extremely quickly. Crucially, particularly for the prop- 
agation of bubbles, the sweep-up time can be significantly 
smaller than a single crossing time; this is relevant because 
a purely hydrodynamic bubble will generally self-disrupt into 
a torus, or smaller fragments in a turbulent medium, in on 
order a crossing time (Robinson et al. 2004; Pavlovski et al. 
2007). 

This buildup of magnetic field will greatly effect the flow 
in the direction of the ambient field lines, and the projectile 
will leave a magnetic bow wave behind it; by analogy with 
other similar bow waves, we expect it to have an opening an- 
gle of tan6' w va/u. In the plane perpendicular to the am- 
bient magnetic field, however, the magnetic field will have a 
much less direct effect as field lines can simply slip around 
the projectile and instabilities can occur In the potential flow 
simulation, the flow smoothly reattaches at the rear of the pro- 
jectile; however, in this case, vortical motions generated at the 
magnetic contact (where the magnetic pressure and magnetic 
tension force is misaligned with the density gradient) and by 
instabilities at the magnetic interface (which are not stabilized 
in this plane) detach the wake from the object, leaving largely 
two-dimensional vortical motions along the field lines in this 
plane. The resulting expected geometry is shown in a cartoon 
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Fig. 4. — We compare the draped magnetic pressure and ram pressure in 
the plane that is parallel (perpendicular) to the initial magnetic field in the 
top (bottom) panels with a logarithmic color scale. In the parallel plane, the 
overpressure of the magnetic draping layer is only partly compensated by 
a deficit in ram pressure and eventually is responsible for decelerating the 
core due to magnetic tension. In the perpendicular plane, the ram pressure in 
the wake of the core attains much higher values and the draping layer closes 
towards the symmetry axis unlike the parallel plane where a nice opening 
cone forms. Shown here is a zoom-in on a small region of our computational 
domain. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 

in Fig. 3 and for our simulation in Fig. 4. This shows that 
the draping layer becomes dynamically important and fills in 
the deficit of ram pressure. The sum of the magnetic and ram 
pressure shows an over-pressure ahead of the core that leads 
to a deceleration of the projectile. 

3. METHODOLOGY 

In order to understand the full non-linear physics of mag- 
netic draping around a dynamically evolving dense projectile 
moving in a magnetized plasma we perform our analysis in 
two steps. First, we analytically study the properties of the 
flow of an ideally conducting plasma with a frozen-in mag- 
netic field around a sphere to explore the characteristics of the 
magnetic field near the surface of the body. To this end, we 
disregard any possible change in the flow pattern by means of 
the back-reaction of the magnetic field. While the derivation 
of this problem can be found in Appendix A, we summarize 
the key results in this section. In the second step, we compare 
this analytical solution to an MHD adaptive mesh refinement 
simulation and explore it quantitatively. 

3.1. Analytical solution 

The potential flow solution for an incompressible flow 
around a spherical body reads as 

,3 X / ^3 



V = Cr 



MCOSfy+t 



MSinf 



(3) 



where R denotes the radius of the sphere and u is the speed of 
the core through the quiescent ambient fluid. Using this solu- 
tion, we solve for the resulting frozen-in magnetic field while 
neglecting its back-reaction onto the flow. For convenience, 
we show here the approximate solution which is valid near the 
sphere. 



2 3s sinO 
B,.= -Bq\ i 

3 \l R l+cose 



(4) 
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Fig. 5. — Diagram showing tlie geometry of the simulations presented here. 
A spherical projectile with a smooth density profile pmax(l +cos(7rr/R))/2 is 
sent in the +z direction with an initial velocity v through an ambient medium 
with density po and a uniform magnetic field pointed in the +y direction. 
The magnetic field strength 'turns on' through the domain with a tanh-profile 
in the direction of motion of the projectile, as indicated by the shading of 
the box; this allows us to start the core in an essentially field-free region and 
smoothly enter the magnetized region. Periodic boundary conditions are used 
in the directions perpendicular' to the direction of motion, and zero-gradient 
'outflow' conditions ai'e used in the z direction. 



(5) 




(6) 



■where ■we introduced a radial coordinate from the surface of 
the sphere, namely s = r-R. These approximate solutions uni- 
formly describe the field near the sphere with respect to the 
angle 9. As described in Appendix A, the energy density of 
the magnetic field forming in the wake behind the body is pre- 
dicted to diverge logarithmically. We point out that the valid- 
ity of the potential flow solution heavily relies on the smooth 
irrotational fluid solution where the magnetic back-reaction is 
negligible. We will see that these assumptions are naturally 
violated in the wake. 

3.2. Numerical solution 
3.2.1. Setup 

The simulations presented in this paper are set up as shown 
in Fig. 5. For clarity of understanding the physical pic- 
ture, we consider only the magnetohydrodynamics (MHD); 
no external- or self-gravity is considered, and we defer other 
physics such as self-consistent inclusion of thermal conduc- 
tivity to future work. In this report, we also consider only the 
magnetic field of the external medium, and assume that it is 
uniform over the scales of interest here. 

In the code units we consider here, the ambient material 
has a density of po = 1, and a gas pressure P = I. The 
(unmagnetized) fiducial projectile has a radius that we vary 
in our runs between R = 0.5 and 2, and a maximum den- 
sity of pmax = 750. With the density profile chosen p{r) = 
Pmax(l +cos(7rr//?))/2, the mean density of the projectile is 
(1 -6/7r^)/9max/2 ~ 0.2pn,ax- Both the ambient and projec- 
tile material are treated as ideal, perfectly conducting fluids 
with ratio of specific heats 7 = 5/3, and so the adiabatic sound 
speed in the ambient medium is y^5/3. The pressure inside 



the projectile is chosen so that the material is initially in pres- 
sure equilibrium. 

The projectile initial velocity is typically chosen to be 1 /4, 
for a Mach number into the the ambient medium of w 0.32. 
The simulation in the transverse directions range from [-4,4], 
and in the direction of motion of the projectile ranges from 
[0,28] for an aspect ratio of 2:7; in most of the simulations 
with projectiles larger than the fiducial R = I, the domain 
size is increased proportionately. The initial magnetic field 
strength can be defined in terms of ao = pu'/Psf,; a typical 
value used in these simulations is 25/4, or Pb„ = 1/100. Pe- 
riodic boundary conditions are used in the directions perpen- 
dicular to the direction of motion, and zero-gradient 'outflow' 
conditions are used in the z direction. Experiments with dif- 
ferent horizontal boundary conditions produced no major dif- 
ferences in results. 

The projectile fluid is initially tagged with a passive scalar, 
so that the material corresponding to the projectile can be 
traced throughout the simulation. 

3.2.2. Code Choice 

As can be seen from analytic arguments, and is shown 
in some detail in §2, two features characterize the problem 
of magnetic draping: the formation of a narrow strongly- 
magnetized layer, and the relative simplicity of the dynam- 
ics, in that a potential flow solution with only magnetic field 
kinematics captures much of the problem, lacking only the 
magnetic field back-reaction. 

Because of the separation of scales (a relatively large object 
moving through an ambient medium and a relatively small 
layer forming around it), the highest resolution requirements 
would impose a large cost on the simulations if the resolution 
had to be everywhere uniform; indeed, it is only a small por- 
tion of the simulation domain which needs to be resolved at 
the highest level. This is especially true since the simulations 
we will need to perform are three-dimensional (as we will 
see in the next sub-section, it is impossible to do meaningful 
simulations of magnetic draping in two dimensions). Thus, a 
simulation code which allows some adaptivity of meshing is 
extremely helpful for approaching this problem. 

The relatively straightforward magnetic field dynamics 
means that, unlike in problems of (for instance) studying the 
details of MHD turbulence, we do not require high-order fi- 
nite difference methods; this is particularly true because of the 
sharpness of the thin layers and the large density gradients in 
this problem. Instead, an MHD solver which can accurately 
deal with sharp gradients is valued. 

As a result of the importance of AMR for these simulations, 
the code we've chosen to perform these simulations with is the 
Flash code (Fryxell et al. 2000; Calder et al. 2002). Flash 
is an adaptive-mesh general purpose astrophysical hydrody- 
namics code which is publicly available'. The MHD solver 
we use here is a dimensionally-split second-order accurate 
8-wave Godunov-type solver which is described in more de- 
tail in Powell et al. (1999). The smallness of spurious mag- 
netic monopoles is ensured by a diffusion-type 'div-B' clean 
operation. This diffusive cleaning approach can be prob- 
lematic near strong shocks, where diffusion cannot operate 
quickly enough; however, no such shocks occur in these simu- 
lations. Flash has been often used for related problems such 
as hot and magnetized bubbles in the intercluster medium 
{e.g., Robinson et al. 2004; Pavlovski et al. 2007; Heinz et al. 

' http://flash.uchicago.edu 
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2006; Roediger et al. 2006; Gardini 2006; Pope et al. 2005; 
Dalla Vecchia et al. 2004; Heinz et al. 2003; Bruggen 2003; 
Briiggen and Kaiser 2002). 

3.2.3. Parameters 

Performing simulations of draping over a projectile with an 
explicit hydrodynamics code (so that compressibility effects 
will be included, for ease of comparison with later, super- 
sonic, work) with a finite resolution places some restrictions 
on the range of parameter space which can be explored. 

For simulating these cases with no leading shock, we re- 
quire that the velocity of the projectile, u, be less than the 
sound speed in the ambient fluid - but to take a reasonable 
number of timesteps (avoiding computational expense and 
spurious diffusion) requires that the projectile velocity remain 
of the order of the sound speed; thus u < Cj, or pou^ < ^P, 
where po and P are the unperturbed density and pressure of 
the ambient medium. For the hydrodynamics of draping to 
be realistic, the magnetic pressure in the fluid must be signifi- 
cantly less than the gas pressure, Pb„ ^ P. Finally, resolution 
requirements for resolving the thickness of the magnetic layer 
will put some constraint on the thickness of the magnetic layer 
I > R/N from Eq. 2, with being the number of points which 
resolve the radius of the projectile; typically the size of the do- 
main (if at full resolution) will be 8N x8Nx 28N. This con- 
straint, expressed in terms of the relevant pressures (ram pres- 
sure and initial magnetic pressure) Pb„ > 3a(pou^)/N. Com- 
bined, these constraints give 



N 



(7) 



For a given u - which is more or less arbitrary, fixed to be 
near the (arbitrary) sound speed - there is thus a relatively 
narrow range of initial magnetic pressures in terms of the ram 
pressure of the ambient material onto the projectile which can 
be efficiently simulated. As we will see, a w 2, and for the 
simulations presented here, ~ 32-64, meaning we are con- 
strained to study roughly that part of parameter space where 

3.2.4. Comparison to Previous Numerical Work 

Having described the setup of the simulations we perform 
for this work, it is worth comparing this to previous numeri- 
cal work. One body of work {e.g., Miniati et al. 1999; Ore- 
gon et al. 1999, 2000) considered a very similar problem in 
the context of the ISM, where a dense cloud may be moving 
through (for instance) the magnetic field associated with the 
Oalactic spiral arms. In this context, a cloud of cold-phase 
ISM gas may very well not be globally self-gravitating (or 
only weakly so) and so a roughly 'top-hat' density function is 
used to describe the cloud. In addition, as appropriate to the 
ISM, supersonic motion was imposed on the ambient mate- 
rial, with the cloud experiencing a high-Mach number shock 
and an ambient magnetic field at the same time; this might 
be appropriate for instance as a cold-phase cloud encounters 
a hot-phase medium for the first time. In our work, we are 
considering subsonic motion initially, so that one significant 
difference is the flow speed of the draping material, in addi- 
tion to the much flatter density profile used in the previous 
work. While we defer the supersonic case to future work, we 
note that in our context, for example in the case of (for in- 
stance) a core falling into a cluster, we would expect the core 
to slowly accelerate before encountering significant magnetic 



field, so the double transient of being shocked and exposed to 
magnetic field at the same time would not be appropriate for 
the case we wish to consider. 

Other more recent work, which in large part motivated the 
research presented here, considers squarely the case we are 
interested in - magnetic draping in cluster environments Asai 
et al. (2004, 2005, 2007). In these works, the authors simi- 
larly consider a spherical projectile moving through an am- 
bient magnetized medium. These simulations include more 
physics (e.g. anisotropic thermal conduction, tangled mag- 
netic fields) than we consider here, as the goal of these works 
was to find out if the draping of the magnetic field could sig- 
nificantly suppress thermal conduction into a cold core. These 
simulations show that magnetic draping is potentially impor- 
tant for suppressing thermal conduction across the interface, 
despite the fact that their draping layers are under-resolved by 
at least a factor of 100 (cf. Eqns. 2 and 7). 

In a slightly different cluster environment, MHD effects can 
also stabilize the morphologies of buoyantly rising radio bub- 
bles from AGN in cluster centers. Ruszkowski et al. (2007b) 
carefully set up magnetically isolated bubbles and tangled 
cluster magnetic fields in the ambient medium as required for 
causally or temporally disconnected generating processes of 
these fields. They showed that magnetic draping of cluster 
magnetic fields possibly supported by helicity of magnetic 
fields internal to the bubble can strongly alter the dynamics 
and suppress hydrodynamic instabilities; thus resembling the 
observed morphologically intact X-ray cavities. 

Knowing that under some circumstances thermal conduc- 
tion can be suppressed by magnetic draping and hydrody- 
namic stability can be enhanced, it becomes worthwhile to un- 
derstand the process of draping itself in some detail, in which 
case the additional micro-physics (at this stage) makes under- 
standing this basic process more, not less, difficult. Thus we 
consider here a case with no explicit micro-physical conduc- 
tion and simplified magnetic field geometry. With this choice, 
our simulations remain scale free and can be applied to var- 
ious astrophysical environments. Additionally, our simula- 
tions can be directly compared to analytics and we are able 
to carefully examine the process of the draping itself. Fur- 
ther, the more simplified physics along with the use of AMR 
allows us to run simulations for significantly longer times as 
the more complicated earlier simulations. 

3.2.5. Two Dimensional Results 

In two dimensions, the imposition of a symmetry greatly 
limits possible magnetic field geometries. In an axisymmetric 
geometry, the only meaningful uniform field geometry is par- 
allel to the axis of symmetry, which in this case would also 
be the direction of motion of the projectile; in this somewhat 
artificial case magnetic field could somewhat constrain a pro- 
jectile (or a bubble; Robinson et al. (2004)) but draping could 
not occur 

In planar symmetry, the field can have components out of 
the plane, in the plane parallel to the direction of motion, or 
in the plane perpendicular to the direction of motion. A com- 
ponent out of the plane will only have the dynamical effect 
of adding to an effective gas pressure {{e.g., Chandrasekhar 
1981)); the component along the direction of motion of the 
projectile cannot be draped. 

Previous work {e.g., Asai et al. 2005) has examined the case 
with two dimensional planar symmetry with a magnetic field 
in the plane of the simulation and perpendicular to the direc- 
tion of motion of the projectile. However, in this case, field 
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Run R 



A 


1 

2 


B 


1 


C 


2 


D 


1 


E 


1 


F 


I 


G 


1 


H 


2 



50 

I 

Too 
1 

lyo 



6.25 


64 


6.25 


64 


6.25 


32 


3.125 


64 


1.5625 


32 


6.25 


32 


25 


32 


15.625 


128 



TABLE 1 

Details of 3-dimensional simulations run for this work. 
Simulations were run with an ambient density and pressure 

of 1 IN CODE units, AND 7 = 5/3. SIMULATIONS WERE RUN UNTIL 
MAXIMUM MAGNETIC FIELD ON STAGNATION LINE WAS 
APPROXIMATELY CONSTANT, TYPICALLY 40-80 TIME UNITS. 



Fig. 6. — The above is a plot of the projectile position (calculated here 
by the maximum height at which there is significant projectile material at any 
given time) over time in a two dimensional draping simulation. At about time 
75, the projectile is actually bounced back under the extremely strong mag- 
netic tension which in two dimensions must grow ahead of the projectile. The 
top panel shows B- in a closeup of the simulation domain at three represen- 
tative times during the bounce, with a white contour indicating the position 
of the original projectile material. This figure has been somewhat degraded. 
For a PDF version of this manuscript with higher resolution figures and inter- 
active 3d graphics, please see http://www.cita.utoronto.ca/~ljdursi/draping/. 

lines cannot cannot slip around the projectile, and so as more 
and more magnetic field gets swept up by the projectile, mag- 
netic tension grows mono tonic ally and linearly ahead of the 
projectile until the forces becomes comparable not only to the 
ram pressure seen by the projectile of the ambient medium, 
but to the ram pressure of the projectile as seen by the ambient 
medium. At this point, the projectile trajectory is reversed. A 
figure describing this is shown in Fig. 6, where the magnetic 
tension forces are seen to compress the projectile (with mean 
density w 150 times that of the magnetized medium) before 
repelling it. 

This outcome is hinted at in Fig. 4 of Asai et al. (2005), 
where in 2d the magnetic energy increases linearly and with- 
out bound, while the 3d models reach a maximum magnetic 
energy. 

3.2.6. Three Dimensional Simulations 

A listing of the eight main runs done for this work are 
shown in Table 1, and the basic setup follows the discussion 
earlier The parameters varied are the size of the projectile, its 
velocity, and the strength of the ambient magnetic field. Other 
runs (the equivalent of run B but with half the resolution, or 
with the same resolution but differing boundary conditions) 
were run to confirm that the results did not change; they are 
not listed here. 

To show that these runs were producing results independent 
of resolution, the maximum magnetic field strength along the 
stagnation line for all the runs with the same Pgo ^nd puf^ are 
plotted versus time in Fig. 7. The maximum field strength is 
sensitive to the resolution of the magnetic field layer, but we 
see here that varying the resolution by a factor of two does not 
effect the results, as the field layer is adequately resolved (but 
only marginally in the case of R/Ax = 32). We also see, as 
we'd expect from the discussions in §2, that the field strength 
in the layer does not depend on the size of the core. 




I ■ ' ■ ■ I ■ ' ■ ■ I ' ' ■ ■ I ■ ■ ■ ■ I ' ' ■ ■ ■ ' ■ ■ ■ ■ ' ■ . 
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Fig. 7. — A plot of the maximum magnetic field strength along the stagna- 
tion line for runs with three different projectile sizes and two different resolu- 
tions across the projectile, all with the same velocity into the ambient medium 
and the same ambient magnetic field. Plotted is the maximum magnitude of 
the magnetic field along the stagnation line vs the time (scaled to the crossing 
time of the projectile). The projectile initially sits in an unmagnetized region. 
The maximum magnetic field strength along the stagnation line is a sensi- 
tive measure of whether the structure of the draped magnetic field is being 
resolved; we see here clear evidence that with the resolution used in this sim- 
ulation the draped layer is being adequately resolved. In the low-resolution 
R = i simulation, which also was run in a somewhat smaller box, towards 
the end of the run the draping layer begins to leave the top of the simulation 
domain, leading to the sudden rapid drop in magnetic field strength. 

The magnetic layer in Run G is under-resolved; while the 
value of R/ Ax is the same as other runs, the velocity is higher, 
so that by Eq. 2 the layer is thinner We include this run be- 
cause it demonstrates certain robustness of results; although 
the layer structure is not adequately resolved at the stagnation 
point, other global properties of the magnetic layer (geometry 
and dynamical effects) otherwise remain robust. 

Indeed, one should be careful about what one means by 're- 
solved' . This discussion should not be taken to mean that the 
other simulations are in all respects resolved. In particular, as 
we will see in §4 and §5, and as suggested by Dursi (2007), 
the flow around the bubble in the xz plane {e.g., transverse to 
the initial magnetic field) is unstable to Kelvin-Helmholtz and 
Rayleigh-Taylor instabilities, not stabilized by the presence 
of magnetic field. Since we have not prescribed any small- 
scale dissipative physics in these simulations, these instabil- 
ities will never properly converge with increased resolution 
(e.g., Calder et al. 2002) as new unstable scales are added. 
This can only be corrected by adding small-scale physics, e.g. 
thermal diffusion; this is left for future work, as the relevant 
microphysics is itself a current research problem (Lyutikov 
2007; Schekochihin et al. 2007). While the rate of develop- 
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ment of the instabilities and their properties is very important 
for long-term mixing of the material of the moving object with 
that of the surrounding medium, we restrict ourselves here to 
studying the development of the magnetic layer and its global 
properties. 

A note here on our AMR is in order In Table 1 , we show 
the finest resolution achieved in the simulation. Flash places 
resolution elements based on a second-derivative refinement 
criterion, in an attempt (appropriate for a second-order ac- 
curate code) to reduce error; see Fryxell et al. (2000) for 
more details. Here we refined based on pressure, density, 
and composition (e.g., fraction of projectile material); early 
experiments with also refining on magnetic field components 
showed little or no difference. The result is that the entire re- 
gion immediately surrounding the projectile is fully refined, 
and the wake region is nearly fully refined. However, because 
no interesting flow features are occurring in most of the do- 
main at any given time, the savings from using AMR can be 
substantial; for run G, for instance, the very simple, laminar 
initial conditions require only 9.7 x 10'' zones; by time 40, 
this has increased to w 8.1 x 10^ zones, and by time 75, an 
approximately-steady 9.1 x 10^ zones are being used; this is 
to be compared to the 469.8 x 10^ zones, or a factor of 51 
greater, that would be required to resolve the entire domain at 
this same finest resolution. 

A final feature worth noting in our runs is that for run E, 
Ml = 1 /2{pu^/PBo) ~ 0.78 < 1; that is, this run has the pro- 
jectile moving sub-Alfvenically, if only marginally. However, 
because of the field geometry, we will see that this makes es- 
sentially no difference for the draping. This is simply because 
the Alfven speed in the direction of motion of the projectile is 
zero - no component of the magnetic field points in that direc- 
tion. While the exact imposition of this condition in our initial 
conditions is somewhat artificial in this case, it is always true 
that it is only the component of the magnetic field that lies 
transverse to the direction of motion that will be draped. 

4. COMPARISON OF THEORY AND SIMULATIONS 
4.1. Magnetic field along the stagnation line 

The first comparison we make is to the one-dimensional 
predictions made along the stagnation line, for instance in 
Lyutikov (2006), where a very specific prediction is made for 
the ramp up, with a particular functional form, of the magnetic 
field strength given in Eq. 1, and it is not necessarily clear that 
such a prediction will hold when the projectile begins to devi- 
ate significantly from spherical. 

To make this comparison, we extract the magnetic field 
strength along the stagnation line for an output, and fit it to 
the equation 

\B\ Bo 1 

' ' (8) 

P Po ' 



icantly larger than the initial radius of the projectile; in this 
example, the expansion is a relatively modest 15%. 



R 

(:-zo) 



where Bq, po are known, and the fit is for the parameters R, 
which would correspond to the radius of the sphere, and zo, 
which would be the position of the centre of the sphere. Re- 
sults for a typical output are shown in Fig. 8. Not only do the 
fits well represent the behaviour magnetic field strength, but 
they also suggest a physical interpretation for the functional 
forms interpretation even when the projectile becomes signif- 
icantly non-spherical; R becomes the radius of curvature of 
the working surface of the projectile at the stagnation line. As 
the projectile becomes more distorted, R can become signif- 




Fig. 8. — Shown is, left, the magnetic field along the stagnation line 
in the simulation ('+') and a fitted theory prediction, with the two fitting 
parameters being the position of the peak and a radius giving a character- 
istic fall off of the field strength. On the right are cut-planes along and 
across the initial magnetic field of the density of the projectile, with a cir- 
cle of radius and position given by the fit to the magnetic field structure, 
left. The radius given by the fit corresponds with the radius of curvature at 
the working surface of the projectile. Results are taken from run B at time 
t = 38.75; results from other simulations and other times give similarly good 
fits. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 



4.2. Comparison of the velocity field 

To compare the potential flow calculations, done in the 
frame of the spherical body, with those of the numerical simu- 
lations, we transform the numerical simulations into the frame 
of the projectile. Because the projectile slows down over time, 
we do not use the initial velocity mq for this transformation, 
but measure the instantaneous mass-weighted velocity of the 
projectile, by making use of the fact that we are tracking the 
fluid that initially resided in the projectile by use of an ad- 
vected passive scalar, a. Thus we measure the instantaneous 
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Fig. 9. — Comparison of the velocity field of the analytical solution in the 
kinematic approximation (top panels) with our numerical simulation (bot- 
tom panels) in the plane of the initial magnetic field. While the velocity 
fields resemble each other very well in the upper half-space, there are distinct 
differences in the lower half-space. These are due to the non-linear back- 
reaction of the dynamically important magnetic field in the draping layer 
on the MHD flow that generates vorticity in the wake of the projectile (cf. 
§5.4). Shown here and in the next figures is a zoom-in on a small region 
of our computational domain that extends up to 112 length units and is four 
times larger in x and y direction. Note that we symmetrized the color map 
of the tiv -component in order not to be dominated by one slightly larger 
eddy. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 
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Fig. 10. — Comparison of the ram pressure in our numerical simula- 
tion (left panel) with the analytical solution in the kinematic approxima- 
tion (right panel). Ahead of the projectile, the ram pressure resembles an 
exact potential flow behavior up to the draping layer which can be seen 
as a black layer around the projectile with a deficit of hydrodynamical 
pressure. Non-linear magnetic back-reaction of the draping layer causes 
the flow to depart from the potential flow solution and to develop vortic- 
ity. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 

velocity of the projectile as 



{pav,} 
(pa) 



(9) 



The top panels of Fig. 9 show the analytical solution of the 
velocity field around the spherical body with radius R in the 
kinematic approximation. For convenience and to simplify 
the comparison to the magnetic field visualization, we show 
the Cartesian components of the velocity field. At infinity, the 
fluid is characterized by a uniform velocity v = -e~u. The 
quadrupolar flow structure results from the fluid decelerat- 
ing towards the stagnation line, the successive acceleration 
around the sphere until 9 = n/2 and mirroring this behavior 
in the lower half-plane by symmetry. The bottom panels of 
Fig. 9 show the numerical solution of the velocity field around 
an initially spherical projectile that deformed in response to 
the non-linear evolution of the magnetized plasma. The white 
line reflects the 0.9 contour of the 'projectile fluid' and cor- 
responds to an iso-density contour The quadrupolar flow 
structure in the upper half -plane resembles nicely the analytic 
potential flow solution. As the flow approaches the projec- 
tile and surrounds it, there are important differences visible. 
In the analytical solution, the flow accelerates for 9 < tt/2 
and decelerates for larger angles 9. In the numerical solu- 
tion, the magnetic draping layer is stationary with respect to 
the projectile. This causes the flow almost comes to rest in 
the magnetic draping layer. The back-reaction of the mag- 
netic draping layer on the flow casts a 'shadow' on the wake 
of the projectile. It prevents the flow to converge towards the 
symmetry axis and suppresses the deceleration of the flow. In- 
stead, vorticity is generated at the draping layer which will be 
studied in detail in §5.4. The comparison of the ram pressure 
in Fig. 10 underpins this argument. 

^ Note that the apparent grid structure seen in the upper part of our simu- 
lated Uj-component is an artifact of our plotting routine as well as small grid 
noise. The interpolation scheme of the plotting routine falsely interpolates a 
smooth velocity gradient with an entire AMR block while it actually drops 
quickly to the velocity value at infinity. 
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Fig. 1 1 . — Comparison of the magnetic field in our numerical simulation 
(bottom panels) with the Taylor expansion of the analytical solution in the 
kinematic approximation that strictly applies only near the sphere (top pan- 
els). We show the Cartesian components (left to right, .v, v, z) of the magnetic 
field in the plane that is parallel to the initial magnetic field. There is a nice 
agreement between both solutions in the upper half-space, while there ai'e 
again distinct differences in the lower half-space. The magnetic shoulders be- 
hind the projectile can be identified that prevents the draping layer from con- 
tracting towards the symmetry axis. In addition, MHD turbulence starts to de- 
velop in the wake of the projectile. This figure has been somewhat degraded. 
For a PDF version of this manuscript with higher resolution figures and inter- 
active 3d graphics, please see http://www.cita.utoronto.ca/~ljdursi/draping/. 
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Fig . 12. — Same as previous figure, but in the the plane that is perpendicular 
to the initial magnetic field. Shown is, left to right, the x,y,z components. 
As expected from our analytic solutions, the draping layer forms by piling 
up magnetic field lines ahead of the projectile. The irregular magnetic field 
in the wake is generated by the vorticity that is absent by definition in our 
potential flow solution. This figure has been somewhat degraded. For a PDF 
version of this manuscript with higher resolution figures and interactive 3d 
graphics, please see http://www.cita.utoronto.ca/~ljdursi/draping/. 



4.3. Comparison of the magnetic field 

It is instructive to compare the analytic solution of the 
frozen-in magnetic field in the kinematic approximation to the 
numerical solution in the planes that are parallel and perpen- 
dicular to the initial magnetic field. We compare the individ- 
ual Cartesian components of the field (Fig. 1 1 and 12) as well 
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Initial magnetic field in the plane: 
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Fig. 13. — Comparison of the magnetic energy density in our numerical 
simulation (left panel) with the analytical solution in the kinematic approxi- 
mation (light panel). The top (bottom) panels show the plane that is parallel 
(perpendicular) to the initial magnetic field. In the analytical solution there is 
a narrow magnetic layer draped around the spherical body, while in our sim- 
ulations the draping layer peels off behind the projectile due to vorticity gen- 
eration. The geometry of the magnetic draping layer in the upper half-plane 
is very similar in both planes suggesting there an approximately spherical 
symmetry. In the wake of the projectile, the draping layer forms a charac- 
teristic opening angle while the field fines can swipe around the projectile in 
the perpendicular plane and the draping layer closes towards the symmetry 
axis. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 

as the magnetic energy density in Fig. 13. Note that we only 
show a Taylor expansion of the highly complex exact solution 
as derived in Appendix A. Strictly, this solution applies only 
near the sphere with an accuracy to ©((r-/?)^/^) as well as 
for flow lines that have small impact parameters initially at 
infinity. Using a different expansion, we verified that the gen- 
eral solution has the appropriate behavior of the homogeneous 
magnetic field at infinity in the upper half-space pointing to- 
wards the positive y-coordinate axis, i.e. rightwards in Fig. 1 1 . 
As expected, the y-component of the magnetic field increases 
as we approach the sphere since the field lines are moving 
closer to each other In the immediate vicinity of the sphere, 
the B field attains a dipolar z-component as the field lines are 
carried around the sphere with the fluid and causes them to 
bend in reaction to the ram pressure of the sphere. As pointed 
out by Bernikov and Semenov (1980) the magnetic lines of 
force that end at the stagnation point are strongly elongated as 
the swipe around the sphere parallel to the line of flow reach- 
ing from the stagnation point into the rear This leads to the 
unphysical increase of the magnetic field as it approaches the 
line of symmetry in the wake and eventually to a logarithmic 
divergence of the magnetic energy density there. 

In the upper half-plane, the analytic solution matches the 
numerical one closely. Interestingly, in the region behind the 
deformed projectile, a magnetic draping cone develops that 
stems from the dynamically important draping layer that has 
swiped around the sphere and advected downstream the pro- 
jectile. In addition, the magnetic pressure in the wake of the 
projectile is also amplified by a moderate factor of roughly 
five (cf. Fig. 13). We will show further down, that this field is 
generated together with vorticity in the draping layer In the 
parallel plane to the initial magnetic field, the magnetic drap- 
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Fig. 14. — A plot for the 3d runs presented here showing the mag- 
netic pressure on the stagnation line once a steady value had been achieved 
for this quantity as a function of the mean ram pressure {p{u)^) as seen 
by the projectile. Omitted is run G, for which the magnetic layer was 
under-resolved and thus the maximum magnetic field strength in the layer 
falls much lower; however, as we will see, even this under-resolution 
does not strongly effect other global properties of the magnetic drape. 
This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http ://w w w.cita.utoronto . ca/~lj dursi/ draping/. 

ing cone causes the stationary flow not to converge towards 
the symmetry axis and protects the region in the wake against 
the increase of the magnetic energy without bounds. The nu- 
merical solution can qualitatively be obtained by remapping 
the analytic solution for 9 > n/l onto the coordinate along 
the magnetic draping cone. In the perpendicular plane to the 
initial magnetic field, there is even better agreement between 
the analytic and the numerical solution. The magnetic field 
in that plane lies primarily in its initial y-direction. This be- 
haviour can easily be understood in terms of the field lines 
sweeping around the sphere in a laminar flow. Numerically, 
we simulate the response of the geometry of the projectile to 
the hydrodynamics. Vortices in the wake deform the projec- 
tile leading to a cap-geometry and a mushroom shape of the 
y-component of the magnetic field. This implies that the flow 
lines detach from the dense material of the projectile gener- 
ating furthermore vorticity and MHD turbulence in the wake. 
The turbulent field mixes the Cartesian components which can 
be nicely seen in the Fig. 12. The magnetic pressure summa- 
rizes our results nicely showing the draping cone in the paral- 
lel plane and the mushroom shaped magnetic energy density 
in the plane perpendicular to that (cf. Fig. 10). Note that we 
choose the same color scale as derived from the simulations 
which leads to a saturation of the magnetic energy density in 
the kinematic approximation at the contact of the spherical 
body and on the axis in the wake. 

5. CHARACTERISTICS OF MAGNETIC DRAPING 

5.1. Field strength in draping layer 

The kinematic solution predicts the magnetic pressure di- 
verges at the stagnation point, which is clearly unphysical. 
From our discussions in §2 and §4, we expect that the mag- 
netic pressure in the draping layer should be on order pu', at 
which point the magnetic back-reaction begins to strongly ef- 
fect the flow; to first order there is no dependence on other 
parameters, such as background magnetic field. One would 
expect, too, from looking at figures such as Fig. 4 that the 
maximum magnetic pressure should exceed the ram pressure 
by some factor, as the magnetic pressure distribution at the 
head of the drape is responsible for redirecting the flow in the 
plane of the draping. 
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Fig. 15. — As in Fig. 1, but for Run E; that is, witli tlie projec- 
tile's velocity reduced by a factor of one-half (so that u = 0.125 in code 
units). A PDF version of this manuscript with an interactive 3d ver- 
sion of this figure, following Barnes and Fluke (2007), is available at 
http://www.cita.utoronto.ca/~ljdursi/draping/. 




Fig. 16. — As in Fig. 1, but for run G; that is, with the pro- 
jectile's velocity increased by a factor of two (so that u = 0.5 in code 
units). A PDF version of this manuscript with an interactive 3d ver- 
sion of this figure, following Barnes and Fluke (2007), is available at 
http://www.cita.utoronto.ca/~ljdursiydraping/. 

We can test this by plotting, for all our runs, the steady 
maximum magnetic pressure at the stagnation line (the field 
quantity that is easiest to consistently characterize) versus the 
mean ram pressure seen by the projectile, p{u)^, where (m) is 
the mean of the projectile velocity (calculated as in by Eqn. 9) 
during the run, and p is the ambient density. The plot is shown 
in Fig. 14 and verifies our expectation. 

5.2. Opening Angle 

The magnetic bow wave behind the projectile is expected 
to propagate transversely away from the projectile at along 
the field lines, and of course to fall behind the projectile at 
velocity u. This suggests a natural opening angle in the plane 
along the magnetic field, tan 9 = va/u. That the direction of 
the scaling is correct can be determined by qualitative inspec- 
tion of a sequence of 3d renderings of simulation outputs as 
the velocity changes; e.g.. Figs. 15,1,16 for m = 1 /8, 1 /4, 1 /2 
and va fixed at 0.1414. 

Although the field lines are stretched during the draping, it 
is the initial Va that is relevant, as the stretching of the field 
lines in the z-direction do not effect the propagation speed in 
the y-direction. For instance, consider a z-velocity shear in 
y, V = (0,0, y/r), with B = (0,Bo,0). The induction equation 
gives us B = W X (v X B) = (0,0, Bq/t), so that the magnetic 
field is only changed in the z-direction; thus ua,, = Va ■ y = 

VAB-y= (\B\/^/Arp)(B-y)/\B\ = By/V?^ = ^Afi 

One can quantify the agreement with this scaling by mea- 
suring the opening angle for the drapes in our simulations. 
The maxima of magnetic field on either side of the stagnation 
line in the y-z plane are found and tabulated along the z di- 
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Fig. 17. — Plot of magnetic energy density in the y,z plane for simulations 
with R = 0.5 and, left to right, u = 0. 125, 0.25, 0.5: shown with black fines are 
the fitted opening angles of the magnetic draping layer, omitting the region 
including the material from the projectile. The fit slopes (e.g., tan 6) are 1.24, 
0.515, 0.261, and those predicted by «ao/(«) are 1.13, 0.566, 0.283; this 
agreement is within 10%. This figure has been somewhat degraded. For a 
PDF version of this manuscript with higher resolution figures and interactive 
3d graphics, please see http://www.cita.utoronto.ca/~ljdursi/draping/. 




Fig. 18. — A plot for the 3d runs presented here showing the tangent of the 
fit opening angles of the drape in the yz plane versus v/i / (w) . Data from all 
runs are shown. 

rection of the simulations, and - omitting the regions above or 
near the projectile itself, and the region below which the drape 
becomes weaker than transient features in the wake - lines are 
fit, and the slope gives the (half-)opening angle. The results 
of the fitting procedure are shown for the same three simula- 
tions in Fig. 17, and a scatter plot for all are runs are given 
in Fig. 18. The scatter for this quantity, and agreement with 
the prediction, is somewhat worse than for the other quantities 
we consider, possibly because the large-scale geometry of the 
draping is more sensitive to the boundaries and the finite size 
of the domain than other, more local, quantities. 

5.3. Deceleration by magnetic tension 

In the scatter plots presented above, we use the mean ve- 
locity (m) of the projectile over time, because there is a mea- 
surable deceleration of the projectile. An example, for run 
F, is shown in Fig. 19. Before the projectile encounters the 
magnetic field at z = 10, hydrodynamic drag - in principle ei- 
ther (numerical) viscous drag or the the drag force caused by 
the creation of a turbulent wake - is all that can play a role, 
and for the simulations presented here, it is the second which 
dominates. The well-known form for the drag on a sphere is 
Fy) = \ /Ipu^ACu, or in terms of a deceleration, 



3 pu^ 
8 {p^)R 



(10) 



where Co is the drag coefficient, experimentally known to be 
between 0.07-0.5, with 0.5 for a turbulent wake, pp is the 
density and A is the cross-sectional area of the projectile in 
the direction of motion. 

However, once the magnetized region is reached and a mag- 
netic layer built up, then another force acts on the projectile 
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Fig. 19. — A plot showing, with squares, the evolution of projectile ve- 
locity (calculated as in by Eqn.9) over time for run F. Note that the projectile 
encounters the magnetic field at time 20 in these units. Plotted as a thin line 
is the best fit deceleration, li = -5.98 X 10"^, for times greater than 20. 
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Fig. 20. — A plot for the 3d runs presented here showing the measured 
deceleration of the projectile versus the functional form we expect it to take, 
proportional to 3/8pu^ /({pp)R), where (pp) is the mean density of the pro- 
jectile (in code units, ss 150). Because the magnetic field strength and curva- 
ture varies over the draped layer, there is an undetermined geometrical factor 
in the magnitude of the deceleration; we find it here to be approximately 1 .87. 

- the magnetic tension from the stretched field-lines. This 
transition can be seen in Fig. 19 for run F; other runs be- 
have similarly. We see that the deceleration caused by the 
magnetic field draping is actually significantly stronger than 
the hydrodynamic draping. This magnetic tension force is 
Fj = /{AttR); we know the magnetic strength in the draping 
layer scales as pu^ (§5.1) and so we can write the deceleration 
as 



pu 



Cc 



(11) 



where Cg is a geometric term taking into account the fact that 
both the magnetic field strength and radius of curvature of the 
field lines vary over the 'cap' of the projectile, and we have 
chosen to normalize Co so that Eqns. 10 and 11 have the 
same numeric prefactor for convenience in comparison. We 
can test this scaling, and at the same time empirically obtain 
Cg, by plotting the decelerations for our different runs, as is 
done in Fig. 20; we find Cg ~ 1 -87. 

It is interesting to note first that the two deceleration terms 
scale in the same way, so that their relative importance re- 
mains constant; and that said ordering is such that the mag- 
netic tension deceleration is always more important, by a fac- 
tor of « 3.7, for the case of highly turbulent (Re w 1000) hy- 
drodynamic drag of Cd = 0.5. In the case of our simulations, 
we do not have the resolution to achieve that highly turbu- 




FlG. 21. — Plot of streamlines over the R = 2 projectile through a 
13 = 100 medium. Streamlines are calculated in the frame of the mean ve- 
locity of the projectile. The streamlines are coloured by the magnitude of 
velocity, and the plane is once again colored by magnetic energy density. 
At this time, no instabilities have developed in the plane perpendicular to 
the ambient magnetic field, so fluid flows smoothly over the projectile in 
this plane; however, fluid traveling close to the other plane experience a 
gain of vorticity. A PDF version of this manuscript with an interactive 3d 
version of this figure, following Barnes and Fluke (2007), is available at 
http://www.cita.utoronto.ca/~ljdursi/draping/. 
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Fig. 22. — The magnitude of vorticity is shown in color in the plane 
across and along to the initial magnetic field (left and right panel). The 
dotted lines represent contours of magnetic energy density. This fig- 
ure has been somewhat degraded. For a PDF version of this manuscript 
with higher resolution figures and interactive 3d graphics, please see 
http://www.cita.utoronto.ca/~ljdursi/draping/. 

lent state. The effective Reynolds number of our simulations 
can be estimated by examining the hydrodynamic drag, for 
example in the first 20 time units of Fig. 20. This does not 
quite give enough data to make a good reading, so we ran four 
simulations with the fiducial parameters {R = \ ,/S.x/R = 32), 
varying m, (0.1 25, 0.25, 0.5, 0.75) and outputting only u(t). An 
excellent fit to the data is provided by Cd ~ 0.77, which corre- 
sponds to (see, e.g.. Fig. 34 in 1V,§45 of Landau and Lifshitz 
(1987)) a Re of just under 200; even in this more viscous case, 
the magnetic draping deceleration exceeds the hydrodynamic 
deceleration by a factor of 2.5. 

5.4. Vorticity generation 

The flow pattern around a moving body looks as follows for 
large Reynolds numbers. The flow is laminar and reflects a 
potential flow solution in almost the entire volume except for 
a narrow boundary layer and the turbulent wake. The char- 
acteristic drag coefficient decreases as the Reynolds number 
increases (Landau and Lifshitz 1987). This turbulent wake 
exerts a drag on the body that decelerates it as described in 
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written in the form 
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Fig. 23. — Plot of the source terms for the specific vorticity injection 
rate as defined in Eqn. 1 3 for the plane transverse/pai'allel to the initial mag- 
netic field (left/right panels). The sohd white contour shows the boundary 
between projectile material and ambiemt material. The upper panels show 
the contribution of the baroclinic-type term where Vp is not aligned with 
the thermal and magnetic pressure force. Due to the large density gradi- 
ent, this term dominates the vorticity injection at the interface between the 
projectile and the ambient medium. The curl of the magnetic tension force 
seems to be the dominant injection mechanism in the wake. Note the differ- 
ent (linear) colour scales for the planes across and along the initial magnetic 
field. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 



§5.3 and shown in Fig. 19. This figure shows two distinctive 
deceleration regimes where only the turbulent drag is present 
in the initial phase, where the magnetic field has not been 
switched on, and a magnetic tension dominated drag phase 
at later times. 

An independent argument is provided by Fig. 21. In the 
pure hydrodynamic case, we do not expect any statistical 
anisotropy of the flow pattern around the moving body. How- 
ever in our MHD flow, there is an unambiguous anisotropy 
visible for the stream lines. In the plane perpendicular to 
the ambient initial magnetic field, the fluid flows smoothly 
over the projectile with only mild perturbations for stream- 
lines near the boundary layer In the plane of the initial mag- 
netic field where the draping cone forms, the stream lines are 
bend towards the turbulent wake and experience the genera- 
tion of vorticity <jJ = V xv. The magnitude of vorticity in our 
simulations is shown in Fig. 22. Vorticity is generated as the 
fluid enters the region in the draping layer where magnetic 
field lines are slipping around the projectile, in particular in 
the plane transverse to the initial magnetic field. The result- 
ing velocity field can not any more be described by the poten- 
tial flow solution which causes the analytical solution to break 
down at the magnetic draping layer and behind the magnetic 
shoulder. The vorticity in the wake suggests the presence of 
MHD turbulence that might be responsible for stretching and 
amplifying the magnetic field furthermore. 

We are interested how exactly the topology of the magnetic 
draping layer can be responsible for generating vorticity into 
an initially vorticity-free flow pattern. The equation of motion 
for an inviscid and magnetized fluid without gravity may be 



dv ^ „^ „ 
at at 

( 1 



(12) 



where we define the convective derivative in the first step 
and applied V x B = 47r/ in the last step. The first term on 
the right-hand side describes the potential force due to the 
sum of the isotropic thermal pressure P and magnetic pres- 
sure Z?^/(87r), while the second term describes the magnetic 
tension force. Applying the curl operator to Eqn. (12) and 
identifying the vorticity u; = V x u, we arrive at the equation 
governing the evolution of vorticity: 
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This equation describes the condition that the vorticity is 
'frozen' in the plasma if the last two terms are negligible.-' 
Vorticity is necessarily generated, if the curl of the force field 
generated by magnetic tension does not vanish (referred to as 
curl-tension term). Another source of vorticity is given by a 
flow where Vp is not aligned with the potential force due to 
thermal or magnetic pressure as well as the magnetic tension 
force (referred to as baroclinic-type term). Figure 23 stud- 
ies the relative importance of both source terms. Due to the 
large density gradient that develops at the interface between 
the projectile and the ambient medium, the baroclinic-type 
term dominates the vorticity injection at this interface; but it is 
unimportant elsewhere. The curl of the magnetic tension force 
seems to be the dominant injection mechanism in the wake, 
and contributes over a broader spatial range in the magnetic 
draping layer. We caution the reader that we cannot quantify 
the level of vorticity injected by means of a turbulent bound- 
ary layer and refer to our phenomenological argument at the 
beginning of this section that clearly indicates the importance 
of the magnetic draping layer for the vorticity injection. 

6. INSTABILITIES 

The magnetic tension force as well as the magnetic layer 
geometry has implications for the instabilities experienced by 
the projectile. In §2 and §5, we saw that the flow in the plane 
parallel to the initial magnetic field is stable and the hydro- 
dynamic instabilities are suppressed by the magnetic draping 
layer (suggested by Dursi 2007). In contrast, the flow around 
the projectile in the plane transverse to the initial magnetic 
field is unstable to Kelvin-Helmholtz and Rayleigh-Taylor in- 
stabilities, not stabilized by the presence of magnetic field. 
We will show that the Kelvin-Helmholtz instability remains 

' This can be seen by considering the evolution of an infinitesimal vector 
5x connecting two neighboring fluid parcels, as the fluid moves with the ve- 
locity field. The point initially at position x at time t will be displaced to the 
position jr + i;(jr)Af at time t + At. The neighboring point initially aXx + Sx 
at time t will be displaced to the position x + 5x + v{x + Sx)At at time t + At. 
Hence this 'frozen' cotmecting line evolves according to 



dt 



(5x) = (5x ■ V)v. 



(14) 



which resembles Eqn. (13) if we neglect the last two terms and identify Sx = 
£U)/p initially, where e > is a small quantity. Since the differential equation 
is true for any time, the same relation will hold for all times for the vorticity. 
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Fig. 24. — Kelvin-Helmholtz instability dissolves the projectile in the 
plane perpendicular to the initially homogeneous magnetic field. The flow 
is accelerated at the density enhancements of the stripped material due 
to the Bernoulli effect. In the wake of the projectile, there is a char- 
acteristic length scale of ~ 10 length units between the striped mate- 
rial which corresponds to an unstable mode with a wavelength of 2/3R. 
This figure has l3een somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 

Stronger and leads to gradual disruption of the projectile, al- 
though the impact of the Rayleigh Taylor instability in our 
MHD case is greater than the purely hydrodynamical case be- 
cause of the greater deceleration. For an homogeneous ini- 
tial magnetic field the induced vorticity remains largely two- 
dimensional. 

The projectile is being decelerated by magnetic tension 
as shown in §5.3. This makes the projectile subject to 
Rayleigh-Taylor instabilities with the growth rate (Chan- 
drasekhar 1981) 

2 {Pp)-Po ■ , 32TrpoCcu^ k SirpoCcu^ 



{Pp)+PO 



(pp) R^ko- 4(/,p) 



where ko = 2n/R defines the the smallest wavenumber of the 
system. We also neglect viscosity and surface tension effects, 
and we work in the limit where (pp) ^ po. 

The flow around the projectile causes a shear at the inter- 
face of the projectile that can get non-linear by means of the 
Kelvin-Helmholtz instability and has the growth rate (Chan- 
drasekhar 1981) 
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k ^3ttu 
(Pp) ^0 ~ R \l {pp 



(16) 

Here we neglect viscosity and self-gravity of the projectile 
and apply the maximal velocity shear from the potential flow 
solution around a spherical body, v = 3/2ue0, which is valid 
atr = /?and6' = 7r/2. 

Which instability will eventually dominate and set the rel- 
evant timescale? It turns out that the ratio of growth rates is 
independent of the projectile properties and only depends on 
the wave number of the considered mode. 



127r k ^ Utt 
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where from the previous section, Cc ~ 1-87 takes into account 
the fact that both the magnetic field strength and radius of 
curvature of the field lines vary over the 'cap' of the projectile. 
The largest length scale of the problem is given by the size of 
the projectile in the direction of motion and sets the largest 
timescale of the problem. 
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Thus, we expect the Kelvin-Helmholtz instability in the plane 
transverse to the initial magnetic field to be responsible for the 
eventual disintegration of the projectile. These considerations 



allow us to estimate the associated time- and length-scale on 
which we expect to see the projectile material in the boundary 
layer to become unstable. 
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(19) 

in terms of the length units in the code. This explains nicely 
the instability features in the wake of Fig. 24 that appear every 
10 length units and indicate that a mode that is slightly smaller 
than the projectile dimension is becoming unstable and leads 
to a deposition of projectile material. 

7. DISCUSSION AND LIMITATIONS 

We have investigated in detail the rapid formation of a mag- 
netic draping layer over a projectile, and examined some of 
the immediate dynamical consequences. It is worth consid- 
ering how well these insights continue to hold over longer 
timescales, and whether the draped field can offer much pro- 
tection over significant distances. 

While details of how mixing might take place will depend 
sensitively on the structure of the object in question, one re- 
quirement for a projectile to mix significantly into the sur- 
rounding medium will be for the projectile to sweep past on 
order its own mass in the ambient medium; only then will 
there have been enough shear to significantly disrupt the mov- 
ing object. This requires the projectile to traverse a distance 
L ^ ((pp) / P{))R- For the runs considered in previous sections, 
modeling this while continuing to resolve the magnetic drap- 
ing layer would require extremely costly simulations, even 
with AMR. 

However, at the cost of complicating direct comparison 
with previous simulations, one can gain some insight into 
what will happen over longer times by considering those re- 
gions of parameter space which make the computation more 
feasible. In particular, for this section we perform an analog 
to run B made with a maximum projectile density reduced by 
a factor of 10, so that (pp) /po ~ 15. With this reduced den- 
sity contrast, mixing happens more easily and the projectile 
sweeps past its own mass in a computationally approachable 
time. Results from this run are shown in Fig 25, at a time 
when the projectile has approximately swept through its own 
mass of ambient medium. 

In this run, the same features are seen as in previous sec- 
tions; the development of the strong narrow magnetic field 
layer, the opening angle ~ va/u, and the large-scale vortic- 
ity oriented primarily along field lines generated in the wake. 
However, over long times the anisotropy imposed by the di- 
rection preferred by the magnetic field, and as suggested in 
Fig. 3, becomes much more pronounced, as the projectile be- 
comes extremely aspherical; it is greatly flattened along the 
direction of the magnetic field lines. 

The long-time distortion of the projectile by the magnetic 
field - which, again, is initially strongly subthermal (/? ~ 100) 
and a factor of 6 less than the ram pressure seen by the pro- 
jectile - is particularly evident when seen compared to the 
results of the same projectile moving in the absence of a mag- 
netic field, as in Fig. 26. In this case, density is plotted in 
two cut planes at the same time for the simulation with and 
without magnetic field. 

Several features are immediately apparent. The first is the 
significantly different density distributions between the two 
scenarios. The densest material is more contained in the case 
with magnetic draping, but in the plane along the magnetic 
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Fig. 25. — As in Fig. 1, but for the (pp) = 15 ran, at a time 
where the projectile has swept past approximately its own mass in ambi- 
ent fluid. A PDF version of this manuscript with an interactive 3d ver- 
sion of this figure, following Barnes and Fluke (2007), is available at 
http://www.cita.utoronto.ca/~ljdursi/draping/. 




Fig. 26. — Plots of density shown for the (pp)/po = 15 run, at a time 
when the projectile has swept past roughly its own mass of ambient medium. 
Shown at top are simulations with no initial magnetic field, and at bot- 
tom which a magnetic field as in Run B, with /3 = 100 and pou^//'fc.o = 
6.25. Panels on the left are along the v direction (e.g., the direction of 
the initial magnetic field lines in the second case) and along the .v direc- 
tion (perpendicular to the magnetic field lines) on the right. These snap- 
shots are taken at the same simulation time for the two rans; the projec- 
tile in the magnetized case lags because of the deceleration demonstrated 
in § 5.3. This figure has been somewhat degraded. For a PDF version of this 
manuscript with higher resolution figures and interactive 3d graphics, please 
see http://www.cita.utoronto.ca/~ljdursi/draping/. 

field, material which is stripped off is more extended, piling 
up along the draped magnetic contact. The plane across the 
initial magnetic field lines is even more interesting; here the 
stripped material is much more contained, even with the pres- 
ence of the Kelvin-Helmholtz instability. Here stripped mate- 
rial stays almost completely within a cylinder of radius R= I, 
the initial radius of the projectile, along the path of the pro- 
jectile. 

Also evident is that, although the two simulations are exam- 
ined at the same time, the projectile without magnetic fields is 
significantly further ahead in the domain than the simulation 
with the magnetic field; this is the result of the deceleration 



demonstrated in § 5.3. 

In this work we have made several simplifying assumptions 
to allow us to begin to understand the process of the draping. 
We have neglected consideration of the interior structure of 
the projectile, by for instance omitting any self gravity which 
would be relevant for a minor merger. This will effect the 
rate of stripping of material off of the core, and thus long- 
term evolution, but is unlikely to directly effect the draping 
process itself. We have also not considered any gradient of 
properties in the medium the projectile moves through; while 
this again would effect long term behaviour, the set-up of the 
magnetic draping layer occurs so quickly that it is unlikely 
that any background quantities would greatly change over the 
small distances involved. 

We have also omitted explicit treatment of dissipative mi- 
crophysics Lyutikov (2007); Schekochihin et al. (2007). Fol- 
lowing up with simulations which included these effects self- 
consistency will be important for examining in detail the re- 
sulting sharpness of the cold fronts (as done, for example, 
by Asai et al. (2007)) and the different dissipation physics 
may also well effect the long time behaviour of mixing. How- 
ever, the initial draping layer is set up, and its properties are 
determined, on timescales much faster than the dissipative 
timescales, so these results will be unaffected. 

We have also considered here only subsonic motions 
through the ambient medium. Many of the astrophysical pro- 
cesses where draping is relevant can be supersonic, and so 
an important next step is to consider this case, where a bow 
shock will occur before the magnetic draping layer. While the 
bow shock will almost certainly be well separated from the 
magnetic draping layer, the shock will affect both the geome- 
try of the flow onto the draping layer and greatly amplify the 
importance of the thermal pressure. On the other hand, be- 
hind the bow shock the flow will be subsonic, so much of the 
discussion here will directly apply. Similar detailed studies 
of draping in the supersonic case will be considered in future 
work. Perhaps more seriously, we have considered here only 
the simplest case of an initial magnetic field uniform over the 
scales of interest. It will be necessary to consider more real- 
istic field geometries. This, too, is being considered in future 
work, and will require much more careful treatment of the 
detailed magnetic structure of the field. 

8. CONCLUSION 

This work aims at understanding the morphology and the 
dynamical properties of magnetic draping to set a solid 
ground for its astrophysical applications. A core, bullet, or 
bubble that moves super-alfvenically in even a very weakly 
magnetized plasma necessarily sweeps up enough magnetic 
field to build up a dynamically important sheath around the 
object; the layer's strength is set by a competition between 
'plowing up' of field and field lines slipping around the core, 
and to first order depends only on the ram pressure seen by the 
moving object. This layer is developed very quickly, poten- 
tially faster than a crossing time of the projectile. The energy 
density in the draped layer, at its maximum, exceeds the ram 
pressure by a factor of two, necessary to anisotropically redi- 
rect the flow. This effect has important implications for galaxy 
cluster physics as it suppresses hydrodynamic instabilities at 
the interface of AGN bubbles. It naturally explains so-called 
'cold fronts' by keeping temperature and density interfaces of 
merging cores sharp that would otherwise be smoothed out 
by thermal conduction and diffusion. Other important astro- 
physical impUcations of this effect include draping of the solar 
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wind magnetic field at the bow shock of the Earth as well as 
pulsar wind nebulae. 

In this paper, we compare a simplified analytical solution 
of the problem that neglects the back-reaction of the dynam- 
ically important magnetic field on the potential flow with a 
high-resolution AMR simulation and find very good agree- 
ment between both solutions in the region ahead of the bubble. 
Non-linear back-reaction of the magnetic field in the drap- 
ing layer necessarily implies the generation of vorticity in the 
flow. The induced vortices in the wake deform the core hydro- 
dynamically and eventually cause the magnetic sheath to peel 
off. There is a strong indication that the vorticity generation 
is responsible for the intermittent amplification and stretch- 
ing of the magnetic fields as well as the injection of MHD 
turbulence in the wake of the core. If this withstands further 
critical analysis, this mechanism might have profound astro- 
physical implications for the amplification and generation of 
large-scale magnetic fields in the inter-galactic medium. 

The magnetic layer, once fully developed, has a character- 
istic geometry which we have shown here to be roughly con- 
ical in the plane along the magnetic field lines, with opening 
angle 9 ~ arctan(i;^/M), and remains contained in the perpen- 
dicular plane, with the Kelvin-Helmholtz instability acting on 
the object interface. The magnetic tension in the layer signif- 
icantly decelerates the object, dominating over any hydrody- 
namic drag. 

Over long times, the anisotropy imposed by the field - de- 
spite the fact that the field is initially highly subthermal and 
with an energy density significantly less than the kinetic en- 
ergy of the ambient medium in the frame as the projectile - 
can significantly distort the projectile, keeping it significantly 
more constrained in the plane perpendicular to the ambient 
field, and keeping any stripped material inside the drape. 

If a magnetic draping layer such as generated in these sim- 
ulations were astrophysically observable, it would be possi- 
ble to get independent measurements of the magnetic field 
strength in the ambient medium provided the local gas den- 
sity and the velocity of the moving object is known. The 
tools are both the opening angle of the drape and the thick- 
ness of the magnetic field layer - but not through the layer's 
field strength. Similarly, for a known magnetic field strength, 
we would have an alternate measure of the velocity of the pro- 
jectile. 
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APPENDIX 



DERIVATION OF THE MHD FLOW AROUND A MOVING BODY 

The exact MHD solution 

The full non-linear solution of the MHD flow around a moving and dynamically evolving body is extremely complex because 
of its significant three-dimensionality and non-linearity. Owing to the range of scales involved this problem is perfectly suited 
for an MHD adaptive mesh refinement simulation. In order to gain credibility and improve our understanding of the properties of 
the numerical solution including its scaling behavior, we solve the problem of an ideally conducting plasma around a spherical 
body analytically. To this end, we solve for the flow of a plasma with a frozen-in magnetic field around a sphere to explore the 
characteristics of the magnetic field near the surface of the body. We disregard any possible change in the flow pattern by means 
of the back-reaction of the magnetic field. The same problem has been investigated by Bernikov and Semenov (1980) who find 
that the energy density of the magnetic field forming in the wake behind the body diverges logarithmically. In passing by we 
correct the misconception that lead to this unphysical behavior of their solution and derive a criterion for the breakdown of our 
simplified analytical solution that we then successfully apply to our numerical solution. 

The governing equations of ideal MHD with infinity conductivity are given by 



curl(t;xB) = and divB = 0. 



(Al) 



We solve this system of equations outside the sphere for a given velocity field that is derived for a viscous and incompressible 
flow around the sphere. Without loss of generality, we choose the origin of our spherical coordinate system to coincide with the 
center of the sphere with radius R (Fig. 1) and the z-axis being anti-parallel to the fluid velocity at infinity. The potential flow 
solution of the pure hydrodynamical problem reads in spherical coordinates as follows (Kotschin et al. 1954): 



v=er\ — - 1 ) ucos9 + ee [ ^ + 1 ) usinO = -u + [3e riu ■ e ,) - u] , 
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where we employed the coordinate independent representation of the homogeneous fluid velocity at infinity in the second step, 
u = Criu ■er)+ee{u -ee)- Since for any stream line holds dr/vr = rd9/ve, we can thus derive the equation of the line of flow using 
Stoke's method of the stream function 



p = rsint 



R^ 
1--, 



(A3) 



where p is the impact parameter of the given line of the flow from the z-coordinate axis on an infinitely distant plane in the left 
half-space. We assume a homogeneous magnetic field at infinity in the left half-space pointing towards the positive y-coordinate 
axis yielding the boundary conditions for B: 



Br\^=Bosmesmc 
Writing Eqns. (Al) for the components yields 
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By substituting (A2) into (A7) we obtain the equation for 
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(A5) 
(A6) 
(A7) 
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(A9) 



where vg/irv,-) = -tan9(2r^ +R^)/[2r{r^ -R^)]. Equation (A9) is a linear inhomogeneous first-order partial differential equation 
which can be solved by the method of characteristics. We take r as parameter in the characteristic equations and express the 
variables 9 and (j> in terms of r, using 



dB^ _ dB^ dB^ 80 dt dB^ 
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Thus, on the line of the flow that is uniquely characterized by its impact parameter p at infinity, we obtain a first order ordinary 
differential equation for B^, 
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Integrating this equation by the separation of variables yields the solution for that contains a constant which is determined 
from the homogeneous magnetic field at infinity (A4), 



B,. 



Bq cos { 
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(A12) 



To determine Br and Bg, we turn to Eqns. (A5) and (A6). By multiplying Eqn. (A5) with r and (A6) with sm9, defining 
K = rsm6{VrBe - vgB,), and combining (A6) and (A5), we obtain the equation for K: 



dr rvr 89 



(A13) 



Equation (A13) can again be solved by the method of characteristics as (A9) yielding K = Kp, where Kp is a constant on each 
flow line that is labeled with its impact parameter p. Determining this constant from Eqn. (A4) and substituting for K and Kp 
their values, we obtain the following equation that relates Br and Bg, 



r sin 9 (VrBg - vgBr) = -puBo sin c 



(A14) 



Substituting Bg, expressed in terms of Br, from Eqn. (A14) and from Eqn. (A12) into the solenoidal condition forB (A8), we 
obtain the equation for Br- Similarly, substituting Br, expressed in terms of Bg and following the same steps, leads to the equation 
for Bg: 
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Both equations can again be solved by the method of characteristics, expressing the variables 9 and <f> in terms of r which we take 
to be the independent parameter along the flow lines and using Eqn. (A3). Note that for a potential flow, the variable is always 
constant on the line of the flow by symmetry. The resulting linear inhomogeneous first-order ordinary differential equations are 
easily solved by an integrating factor that is derived from the homogeneous equations, leading to the solutions for Br and Bg, 



Br 



-R' 



■cost 



2r^+R 



Ci T^osint 



C2 ± 2Bq sin (; 



p{r,9)r"^Ar' 



-R^- p(r, 9Yr'f'^ Vr'^-R^ 



r'\r'^ + 2R^Wr'^-R^dr' 



5 (2r'3 + /;3)2 ^^/3 _ ^3 _ ^(^^ 0)2 



(A17) 



(A18) 



where Ci and C2 are integration constants and ^ is the initial value for which Br and Bg are known. The upper signs refer to the 
region <9 < tt/2, and the lower signs to tt/2 < 9 <tt. 

The approximate MHD solution near the sphere 

We aim at understanding the behavior of the magnetic field in the region near the sphere. To this end, we investigate the 
behavior of the integrals in (A17) and (A18) for small impact parameters and keep only the main terms with respect to p. We 
find that the integral in (A17) diverges at the lower limit logarithmically for n/2 and the integral in (A18) has a linear divergence 
at the lower limit. Thus we will use (A17) in the region Q < 9 <ii /2 and (A18) in the region for 7r/2 < 6* < tt. 

We divide the region of integration into two: the first from 00 to ri where r\ > R is the radius of the sphere on which the 
asymptotic form of the magnetic field changes, and the second from rj to ro, where ro is the radial value of the flow of line under 
consideration for 9 = tt/2. This implies that the following expansions only apply to small impact parameters p with r^ < r^. By 
expanding the integrand of (A17) in powers of 1/r for r > ri > R and in the region < 9 < tt/2, we determine Ci = and we 
recover the homogeneous field at infinity with an accuracy to 0(1/ r). Near the surface ri > r > ro we perform a change of the 
variable to s = r-R. We define si = ri -R and so = ro-R and s varies within io < * < ^'i. The equation of the line of flow (A3) 
has the form p = ^/JsR sm9 with an accuracy to 0(s^^^) and from this we obtain sq = p^ /(3R) for 9 = tt/2 and s = sq. 

The value of B,. in this region will be composed of two terms: the value of the integral in (A17) from 00 to ri with a factor to 
leading order (x s^^^, and the value of the integral from 5i to s, which behaves like 0{s^^^). Neglecting the first term in comparison 
with the second, we obtain for B,- with an accuracy to 0(s^^') or O(p^): 



Br = - 



3s 
~R 



Bo P sin cos 9 



sds 



9(s^-ssQ)y^' 



For impact parameters p with sq < si, we obtain with an accuracy to Oiso/si): 

D 2 /37 sin6' . 

Br = T^o \/ — 7- ^ sm (/). 

3 \ R l+cost* 



(A19) 



(A20) 
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Using Eqn. (A14) leads to the component fie- Thus, Bg and near the sphere are determined by the formulae 



I = Bq sin 



Ba,=B[)Cos 




(A21) 



(A22) 



It turns out that these formulae are also correct for the region 7r/2 < 6 < n as follows from Eqn. (A18). The integral in this 
expression is regular for so, and by computing Bg in the approximation s <^ R for s = so, we find C2 = sin0/?Bo/3. Then Bg 
is equal to (A21) with an accuracy to terms of order 0{s^^^). We obtain (A22) by using (A14). Thus, Eqns. (A20) to (A22) 
uniformly describe the field near the sphere with respect to the angle 9. 

In order to facilitate comparison to our numerical solution, we transform the approximate solution for B given by the compo- 
nents in the spherical coordinate system (A20) to (A22) into Cartesian system yielding 



Bt = Bocos(f)sin(f){l-cos9) \ — f 1 

\ 3s \ R 



By = BQy 



sin- </.(!- cos 6*) ( ^-1 ) + l 



= Bx tan + Bo A 



/R f 2s cos6 
3^ V R l+cost' 



(A23) 
(A24) 
(A25) 



Note that we introduced the radial coordinate from the surface of the sphere s = r-R and that this solution applies only near the 
sphere with an accuracy to 0(s^^') as well as for small impact parameters p with an accuracy to O(so/si). 

Using the method of regularization of the integral in (A17) with respect to the lower limit 6 = tt/2, Bernikov and Semenov 
(1980) investigate the behavior of the magnetic field in the wake of the sphere. They find that, when neglecting a term that scales 
as 0(1/ r), Br is given by 

b''-^, (A26) 
3 P 

which, with proximity to the axis of the wake p ^ 0, leads to an unlimited increase of Br ^ 00. The magnetic lines of force 
that end at the stagnation point are strongly elongated as the swipe around the sphere parallel to the line of flow reaching from 
the stagnation point into the rear. This leads to the unphysical increase of the magnetic field as it approaches the symmetry line. 
While this might be the mathematically correct solution, it leads to a logarithmic divergence of the energy density of the magnetic 
field in the volume near the wall. 
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